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ABSTRACT 

(N ; 

^! As part of the In-Orbit Checkout activities for the Spitzer Space Telescope, 

the IRAC team carried out a deep observation (average integration time ~ 8 
| ! hours) of a field surrounding the bright QSO HS 1700+6416. This field contains 

several hundred z ~ 3 Lyman-break galaxy candidates, and we report here on 
\ their mid-infrared properties, including the IRAC detection rate, flux densities 

and colors, and the results of fitting population synthesis models to the optical, 
{Nj \ near-infrared, and IRAC magnitudes. The results of the model-fitting show that 

^ \ previous optical/near- infrared studies of LBGs were not missing large, hidden old 

stellar populations. The LBG candidates' properties are consistent with those 
of massive, star-forming galaxies at z ~ 3. Other IRAC sources in the same 
\ field have similar properties, so IRAC selection may prove a promising method 

of finding additional high-redshift galaxies. 
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1. Introduction 

Observations of galaxies at high redshift provide important tests of theories of galaxy 
and structure formation: galaxy masses, star formation rates, and the distribution of these 
quantities are all vital information. Some of the best-known high-redshift galaxies are the 
z ~ 3 galaxies discovered with the Lyman-break technique (Steidel & Hamilton 1993). These 
galaxies are forming stars rapidly and are generally thought to be massive enough to be the 
ancestors of today's large galaxies (Adelberger et al. 1998; Steidel et al. 1998; Weatherley 
& Warren 2003). Rest-frame near- infrared observations primarily sample the old stellar 
population in galaxies and are not severely affected by dust extinction. They should provide 
a less biased estimator of the galaxy mass than optical measurements, which are sensitive to 
both internal extinction and emission from young stars which make up little of the galaxy 
mass (Brinchmann & Ellis 2000; Kaufmann & Chariot 1998). The four bands of the Infrared 
Array Camera (IRAC; Fazio et al. 2004a) on the Spitzer Space Telescope (Werner et al. 
2004) probe the rest-frame near-IR at redshifts z = 2 — 4. 

This paper reports Spitzer/IKAC 3-10 /im observations of a Lyman-break galaxy field. 
The goals of the observations were to confirm that IRAC could detect z ~ 3 galaxies and to 
characterize the noise performance in the approach to the confusion limit. Galaxy number 
counts in this field are addressed in a separate paper (Fazio et al. 2004b). IRAC properties 
of z ~ 2 galaxies in this field are also addressed in a separate paper (Steidel et al. 2004, in 
preparation). Throughout this paper we assume a cosmology with H = 70 km s _1 Mpc -1 , 
Vt m = 0.3, and fi A = 0.7. 

2. Observations and data reduction 

The field observed was a Lyman-break galaxy field surrounding the z = 2.7 QSO 
HS 1700+6416, from the survey of Steidel et al. (2004). This particular field was chosen 
because of its existing very deep optical photometry, which allowed selection of high-redshift 
galaxies, and its position within the Spitzer Continuous Viewing Zone. The field also in- 
cludes two galaxy clusters (Abell 2246 at z = 0.25 and an anonymous cluster at z — 0.44), 
whose X-ray emission was discovered during a study of the QSO (Reimers et al. 1997). The 
IRAC observations were centered at coordinates 17 h 01 m 17 s , 64 d 09 m 02 s (J2000), about 3.7' 
southeast of the QSO. The central IRAC pointing was offset in order to slightly lower the 
observed galaxy density and reduce the effects of confusion. The two IRAC fields of view 
see different regions of the sky, so the final 3.6 and 5.8 /im images extend further to the 
southeast while the 4.5 and 8.0 /im images extend to the northwest. The IRAC images were 
taken on (UT) 14-16 October 2003 as part of Spitzer Program ID 620. We used 200-second 
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frametimes in the 3.6-5.8 channels; this observing mode uses 50-second frametimes in 
the 8.0 /iiii band because of the higher sky background. We obtained some additional data 
with 100-second frametimes, but only used the 8.0 /im-band data (which again uses 50 s 
frametimes) to avoid systematics in combining data taken with different frametimes. 

The individual images were processed with the Spitzer Science Center (SSC) pipeline, 
then had their astrometric positions refined with 2MASS sources. The zero levels of the 
5.8/im-band images were adjusted to minimize differences in the overlap areas using the 
overlap consistency option in the SSC processing tools. The images were combined into 
mosaics with our IDL-based code which re-projected the images to correct for distortion and 
combined them with a sigma-clipping algorithm. The resulting mosaics are in the original 
array orientation, with a pixel scale half as large as the original. 1 

Figure 1 shows portions of the central regions of the mosaic images in the 4 IRAC 
bands. Because of dithering, the effective exposure time across the mosaics varies, ranging 
from about 11 hours in the deepest areas to 20 minutes near the edges. The average exposure 
time over the central 5 x 10 arcmin area is 7.8 hours. The detection limits vary with the 
effective exposure time; approximate 5a limits in the central regions of the images are 0.45, 
0.45, 0.8, and 0.9 /iJy (AB magnitudes of 25.0, 25.0, 24.2, 24.1) in the four bands. 

The optical data used for LBG selection are described by Steidel et al. (2004). Although 
these data were obtained for a survey of z = 2 galaxies, the imaging methods are the same 
for both z = 2 and z = 3 selection (only the color criteria differ). Briefly, the optical 
imaging observations were obtained at the William Herschel 4.2m telescope in 2001 May 
and supplemented with [/-band images from Keck I/LRIS-B later that year for the central 
region near the QSO. The la surface brightness limits in the 3 bands (in a 1" aperture) 
are U n (AB) = 28.9, G(AB) = 29.0, and R S (AB) = 28.2. A deep fT s -band image with 
Kab — 23.9(5er) covering the central 9x9' was obtained at the Palomar Hale Telescope 
in 2003 June. Optical spectroscopy of an LRIS-B mask targeting 22 z ~ 3 candidates was 
obtained in 2003 September. Reduction followed the method described by Steidel et al. 
(2003), and 17 redshifts were successfully measured, with z = 2.911. The spectroscopically- 
confirmed galaxies had a median R S (AB) = 24.6. 



The mosaics will be available via http://cfa-www.harvard.edu/irac/publications. 
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3. Analysis 

3.1. Object detection and photometry 

Lyman-break galaxy selection was done using the U n GlZ technique described by Steidel 
et al. (2003). A total of 445 LBG candidates in the 15' x 15' area centered on the QSO 
were selected, to a limiting total magnitude of R S (AB) = 25.5. These objects are typically 
> 10(7 detections in 7Z, 15 — 20a" in G, and very faint or undetected in U n . In general, 
spectroscopic observations of LBG candidates show that about 95% are likely to be bona 
fide z ~ 3 galaxies (Steidel et al. 2003). There are about 110 LBG candidates in the 3.6 
/xm/5.8/nn field of view, about 175 candidates in the 4.5/zm/8.0/zm field of view, and about 
90 in the overlap area. Parts of the final 3.6 and 5.8 /im images extend outside the U n G1Z 
detection images. 

To detect LBG candidates on the IRAC images, we used the image world coordinate 
system to project the candidates' positions into the image plane and performed aperture 
photometry at those positions. We experimented with allowing the photometry algorithm 
to re-center the objects; if the objects moved more than l'/2 in either x or y (indicating 
that the identification was in danger of moving to a nearby object), we used the un-centered 
position. This is important because the IRAC images are crowded, particularly at the shorter 
wavelengths (the IRAC point spread function FWHM is about l'.'S-^.'O). As an indicator 
of detection significance, we used the IRAF/APPHOT photometric uncertainty inside a l'/5 
radius aperture, with sky counts determined in 15-20" annuli. Objects with magnitude 
uncertainty < 0.217 mag (S/N > 5) were considered detected (see Figure 1 for examples). 
The final IRAC detection numbers for the 4 bands are given in Table 1. Most of the LBG 
candidates were detected at 3.6 and 4.5 /xm; the detection rate decreases with increasing 
wavelength to 45% at 8.0/im. Photometry of the LBG candidates was completed by applying 
aperture corrections from the l'/5 measurement aperture to the 12'.'2 aperture used for flux 
calibration. The corrections, derived from the measured IRAC PSF (S.T. Megeath, priv. 
comm.), are —0.52, —0.55, —0.74, and —0.85 mag. 

To generate a sample of IRAC-selected non-LBG objects for comparison to the LBG 
candidates, we used DAOFIND to find objects to a threshold of 5<r above local background. 
The detection images were the IRAC mosaics multiplied by the square root of the coverage 
maps (to give uniform noise across the images) with areas near the image edges, near bright 
stars, or near instrumental artifacts masked out. Detection was performed in the 3.6 and 
4.5 fim band images only, with object positions transferred to the frame of the 5.8 and 
8.0 /im band images for photometry. For a more reliable sample, we included only objects 
detected in both 3.6 and 4.5 /im bands, which restricts the area used to about 42 arcmin 2 . 
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About half of the LBG candidates in the field were detected as part of this sample (the rest 
were too close to brighter objects or masked areas), and were removed. The final non-LBG 
sample contained 1931 objects. Photometry of the non-LBG sample was done in an identical 
manner to that of the LBG candidates. No attempt was made to separate stars and galaxies 
in the non-LBG sample on the basis of spatial extent; the size of the IRAC PSF makes this 
impractical for distant galaxies. The galactic coordinates of the field (/ = 94.3, b = +36.2) 
are such that we do not expect stellar contamination to be a serious problem for faint objects. 

3.2. LBG properties: flux and color distributions 

To compare the optical properties of IRAC-detected and non-detected LBGs we used 
the optical catalog generated for LBG selection. There was no significant difference (K-S 
test probability > 5%) in 7Z between objects detected and not detected in the 3.6 /xm band. 
However, there are so few non-detections in this band (only 10) that this comparison is not 
particularly meaningful. A K-S test of the 1Z magnitude distributions of IRAC detections 
and non-detections in the 4.5, 5.8, and 8.0 jum bands showed that the differences in mean 
magnitudes are statistically significant. The IRAC-detected objects are roughly 0.2-0.4 mag 
brighter in 1Z than the non-detected objects. There is, however, no clear correlation between 
TZ and IRAC magnitudes for the IRAC-detected objects. 

We can also compare the IRAC flux density distributions of the LBG candidates and 
the non-LBG sample. The LBGs have median IRAC magnitudes 0.5-0.7 fainter than the 
non-LBG sample. The spectroscopically-confirmed galaxies are slightly brighter than the 
full candidate list in the 3.6 and 4.5 /im bands and slightly fainter in the longer- wavelength 
bands (but their numbers are small: only six such galaxies have four IRAC-band detections). 
The 33 LBG candidates detected in all four IRAC bands have median AB magnitudes in the 
IRAC bands of 22.7, 22.7, 22.6 and 22.5 (flux densities in /iJy: 3.1, 3.1, 3.3, 3.8). Two-color 
plots are shown in Figure 2. The LBGs are redder than the comparison objects in [3.6]-[4.5], 
for which the color limit is constant with magnitude. Selection effects bias the observed color 
distributions in [3.6]-[5.8] and [3.6]-[8.0]: because of the different sensitivies in the IRAC 
bands, bluer objects are observable only at bright [3.6] magnitudes. 

Figure 2 also shows predicted colors for template SEDs. The templates combine optical 
and near-infrared data from the hyperz package (Bolzonella, Miralles, & Pello 2000) and 
mid-infrared data from Lu et al. (2003). As Papovich, Dickinson & Ferguson (2001) report 
for optical/near- infrared data, the 'irregular galaxy' template SED is the closest match to the 
LBG candidates' colors, particularly the spectroscopically-confirmed objects. Besides being 
known to be at the correct redshift, the spectroscopic galaxies also tend to be brighter, so 
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their colors should be better determined. A substantial number of LBG candidates fall well 
away from the template colors, which could happen for several reasons: photometric error 
(for example, due to contamination by nearby objects), incomplete or inapplicable templates. 

Some objects in the comparison sample have the same fluxes and colors as the LBG can- 
didates, and therefore could be high-redshift galaxies. Using the top two panels of Figure 2, 
we suggest selection criteria for z > 2 galaxies based on both the LBG candidate observations 
and the template tracks. The criteria [3.6]ab — [4-5] ab > —0.15, 1.1 > [3.6]ab — [5.8]ab > 
—0.4, 1.1 > [3.6]ab — [8.0]ab > —0.5, and [3.6]ab > 22 give good separation between low 
and high-redshift galaxy templates, while including about two-thirds of the LBG candidates 
and all but one spectroscopically-confirmed LBG. This color criterion would also pick up 
some M82-like SEDs at 1 < z < 2; however, the bottom panel of the Figure shows that 
adding the criterion R AB — [3.6] < 4.1 + 2.4([3.6]ab — [4-5]ab) should reduce the number 
of such objects. These criteria define a subsample of 86 objects from the non-LBG sample, 
giving an areal density of approximately 2 arcmin" 2 . This is comparable to the number 
density of LBGs in the same area (2.0 arcmin^ 2 ). While the LBG technique is designed to 
select galaxies over a fairly narrow redshift range, the use of IRAC/optical colors may permit 
selection of high-z galaxies without biases toward UV-bright objects, albeit with an as yet 
undetermined contamination rate (see also Huang et al. 2004). Confirmation of such color 
selection with spectroscopic redshifts is of course necessary, and this will be the subject of 
future investigations. 



3.3. Spectral energy distributions and model fitting 

At high redshifts, the IRAC bands sample the rest-frame near-IR light of galaxies. This 
light is expected to be a good tracer of galaxy stellar mass since it comes mainly from low- 
mass older stars. To test this assumption, we compared stellar population synthesis to the 
combination of optical, near-infrared, and IRAC data. Following Shapley et al. (2001), we 
used the population synthesis models of Bruzual & Chariot (2003) to generate synthetic spec- 
tral energy distributions for composite stellar populations with solar metallicity, a Salpeter 
(1955) initial mass function, and star formation rates (SFRs) 

q?(t) = ^ e~* /TSFR (1) 

with tsfr = 10, 50, 100 Myr, and galaxy ages at z = 3 of 10, 50, 100, 200, 500, 1000, 
1200, 1500 and 2000 Myr, corresponding to formation redshifts 3.21 < Zf < 30. As a 
matter of convenience, we used the Bruzual & Chariot (2003) code's internal prescription 
for extinction, derived from Chariot & Fall (2000). Models with four different values of the 
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extinction parameter ry (0.5, 0.75, 1.0, 1.5) were generated. Predictions for observed fluxes 
were made from the model SEDs using the program cm_evolution included in Bruzual & 
Chariot's software distribution. We expect the rest-frame UV photometry to provide a more 
sensitive estimate of the instantaneous SFR in these models, while the rest-frame near-IR 
photometry from IRAC will provide the constraints on the integrated stellar mass. 

The best-fitting models were determined by minimizing x 2 , the sum of the difference 
between model and observations. The Bruzual & Chariot (2003) models have a tabulated 
stellar mass for each age, so the overall best-fit normalization of the model apparent magni- 
tudes fixes the model stellar mass. We fit models to the GRK S +IKAC magnitudes of each of 
the six spectroscopically-confirmed galaxies individually and to a 'composite' SED made from 
the median fluxes of the the 26 unconfirmed candidates (assumed to be at z = 3.00). The 
model SEDs and observations, with model parameters, are shown in Figure 3. Current star 
formation rates ^(t bs) are derived using Equation 1 and requiring J * obs ty(t)dt = (1 — R)M* 
(R = 0.28 is the fraction of formed stellar mass returned to the ISM; Cole et al. 2000) 

Most of the observations are best fit by models with star formation t S fr = 50 — 100 Myr, 
relatively young ages (100-300 Myr), and stellar masses in the range 1.5 — 4 x 10 10 M Q . The 
resulting star formation rates are 7— 33M Q yr _1 . The largest SFR is for D17, a galaxy without 
a measured i^-band magnitude which also has a large best-fit extinction. This object has 
bluer IRAC colors than the predicted SED of the prototype dusty starburst M82 but could 
be intermediate in type between M82 and a nearby irregular. Previous results on SED 
fitting of LBGs with optical and near-infrared observations are presented by Shapley et al. 
(2001) and Papovich, Dickinson & Ferguson (2001); the latter used a much more extensive 
grid of population synthesis models. Our model parameters are well within the range of 
best-fit parameters found by these two groups, indicating that IRAC-detected LBGs are not 
particularly unusual members of the LBG population, and the previous optical/near-infrared 
studies were not missing large, hidden old stellar populations. More sophisticated model- 
fitting, using a wider range of stellar population parameters and extinction prescriptions, 
is clearly possible but beyond the scope of this paper. Such analysis will be pursued in 
subsequent work, using these data and a larger sample of LBGs in the Extended Groth Strip 
survey area. 



4. Summary 

Deep observations with IRAC detect nearly all Rab < 25.5 Lyman-break galaxy can- 
didates at 3.6 and 4.5 /im, and about half of the candidates at 5.8 and 8.0 /im. The color 
and magnitude distributions of LBGs overlap with those of a sample of comparison objects 
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in the IRAC field. This suggests that IRAC selection may allow detection of different types 
of high-redshift galaxies than are found by Lyman-break selection. Spectral energy distribu- 
tion fitting implies that the IRAC-detected LBGs are massive stellar systems with relatively 
recent star formation. The sensitivity of IRAC to old stellar light at these high redshifts 
shows its promise for future work in galaxy formation and evolution. 

We thank M. Ashby for a careful reading of the manuscript, and the referee for helpful 
comments. This work is based on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through Contract Number 
1256790 issued by JPL/Caltech. Support for the IRAC instrument was provided by NASA 
through Contract Number 960541 issued by JPL. M.A.P. acknowledges NASA/LTSA grant 
# NAG5-10777. 

Facilities: Spitzer(IRAC) 

The data from this paper, ads/sa.spitzer^0007127552, ads/sa.spitzer#0007127808, 
ads/sa.spitzer#0007128064, ads/sa.spitzer#0007128320, ads/sa.spitzer#0007128576, 
ads/sa.spitzer#0007475968, ads/sa.spitzer#0007476224, ads/sa.spitzer#0007476480 
is available through the electornic edition. 
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Fig. 1. — IRAC images of the Q1700+64 field in the 4 bands (wavelength labels in the upper 
right of the images). These images show a portion (3.5' x 3.5') of the area covered by all four 
IRAC bands; north is to the upper right and east to the upper left. White circles indicate 
LBGs detected in all four IRAC bands; black circles are LBGs detected in three bands; white 
squares are LBGs detected in two bands. 
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Fig. 2. — Two-color plots of LBG candidates (asterisks), confirmed z ~ 3 LBGs (filled 
squares) and non-LBG sample (dots). Lines are template SEDs for 0.5 < z < 4, with 
redshifts 0.5,1,2, and 4 marked. The different templates are: Sbc spiral (thick solid), irregular 
(dashed), M82 (dotted). Thin solid lines indicate color-selection criteria for high-redshift 
galaxies. 
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Fig. 3. — Spectral energy distributions for spectroscopically-confirmed LBGs and a com- 
posite of the photometric candidates. Low-resolution SEDs for the best-fitting Bruzual 
& Chariot (2003) models are shown; models are normalized to the data at A G bs= 3.6/iin 
(A rest ~ 0.9/im). From shortest to longest wavelength, the photometric points marked are in 
the G, 1Z, K S) 3.6, 4.5, 5.8, and 8.0 fj,m bands. Ages and tsfr are in Myr; stellar masses are 
in M Q . 
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Table 1. IRAC Detections of Lyman break galaxies 



Photometric candidates Confirmed z ~ 3 



IRAC band N(U n GTZ) a N(IRAC) b N(LRIS) a N(IRAC) b 



3.6 /xm 


114 


104 


14 


14 


4.5 /xm 


180 


150 


10 


10 


5.8 /xm 


112 


68 


16 


14 


8.0 /xm 


170 


77 


10 


7 



a Number in field. 

b Number detected with IRAC. 



